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ABSTRACT 

After a short review of experimental foundations of met- 
ric theories of gravity, the choice of general relativity as 
a theory to be used for the routine modeling of Gaia ob- 
servations is justified. General principles of relativistic 
modeling of astronomical observations are then sketched 
and compared to the corresponding Newtonian princi- 
ples. The fundamental reference system - Barycentric 
Celestial Reference System, which has been chosen to 
be the relativistic reference system underlying the future 
Gaia reference frame is presented. Principal relativistic 
effects in each constituent of a relativistic model of astro- 
nomical observations are briefly elucidated. The structure 
of a relativistic model of positional observations which 
can be used as a standard relativistic model for Gaia is 
sketched. The physical meaning of the Gaia reference 
frame is discussed. It is discussed also how Gaia obser- 
vations can be used to verify general relativity. 

Key words: relativity, reference systems, Gaia reference 
frame. 



1. WHY RELATIVITY? 

Reduction scheme of positional observations in Newto- 
nian physics is rather simple. Absolute Euclidean space 
and absolute time of Newtonian physics lead to the exis- 
tence of global preferred coordinates: inertial coordinates 
which are unique up to a constant shift of the origin of the 
time coordinate, constant rotation of spatial axes and a 
shift of the origin of spatial coordinates which is at most 
linear in time. Although already in Newtonian physics 
one can introduce arbitrary coordinates (e.g., some curvi- 
hnear coordinates), the inertial coordinates are certainly 
preferred since the laws of physics look especially simple 
when expressed in an inertial reference system. More- 
over, observed quantities (distances, directions, etc.) are 
directly related to those global inertial coordinates. 

Let us briefly consider the Newtonian scheme of reduc- 
tion of astronomical observations. Figure ^ sketches the 
four constituents of an astronomical observation from the 
point of view of Newtonian physics: (1) motion of the 
observed object, (2) motion of the observer, (3) propaga- 
tion of an electromagnetic signal from the object to the 
observer, and (4) the process of observation. The last two 
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Figure 1. Four parts of an astronomical event from the 
point of view of Newtonian physics: 1) motion of the ob- 
served object; 2) motion of the observer; 3) trajectory of 
an electromagnetic signal from the observed object to the 
observer which is tacitly assumed to be a straight line in 
Newtonian astronomy; 4) the process of observation re- 
sponsible for Newtonian aberration. The coordinate grid 
in the background symbolizes a global inertial reference 
system. 



parts can be formulated in a quite simple way in Newto- 
nian physics. It is normally tacitly assumed here that the 
light rays are straight lines in some inertial coordinates. 
As for "the process of observation", it is responsible for 
the appearance of Newtonian aberration which reflects 
the difference in observed directions to the source by a 
moving observer and by an observer at rest relative to the 
chosen coordinates. 

The goal of Newtonian reduction of astronomical obser- 
vations is to model (to predict) the results of observations 
performed by a fictitious observer (normally situated at 
the origin of the chosen reference system, e.g. at the 
barycenter of the solar system) at some given moment 
of time. One attempts here to correct for all the effects in 
observations which are produced by the motion and the 
position of the real observer (aberration and, e.g., paral- 
lax, respectively) and by the motion of the object (proper 
motion and, possibly, light travel time effects). The struc- 
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ture of a Newtonian reduction scheme does not depend 
on the goal accuracy of reduction and can be described 
as follows: (1) aberration, (2) parallax, (3) proper motion 
and/or light travel time effects. For low accuracies when 
only linear effects from aberration, parallax and proper 
motion are of interest, one could apply the correspond- 
ing corrections in arbitrary order On the contrary, for 
higher accuracies the order of these reductions is impor- 
tant. All parameters of the model, i.e. the coordinates of 
the observer and the object as function of time, are de- 
fined in the chosen inertial reference system. That is, the 
five standard astro me trie parameters of the object (right 
ascension a, declination S, parallax tt, proper motion in 
right ascension fia and proper motion in declination /i^) 
are also defined in the chosen reference system. 

Rapid increase of observational accuracy of astronomi- 
cal observations has already made indispensable to use 
general relativity for modeling of the observational data. 
For many kinds of observations the Newtonian scheme 
sketched above fails to describe observational data with 
the required accuracy. In many cases the deviations from 
the model are several orders of magnitude larger than 
the accuracy of observations. Examples are astrometric 
(geodetic) VLBI observations, lunar laser ranging, radar 
ranging to the planets, experiments with high accuracy 
clocks, GPS observations. It is also widely known and 
accepted that the deviations can be eliminated by using 
Einstein's general theory of relativity (instead of Newto- 
nian physics) for the modeling of observations. 

The accuracy of positional observations to be produced 
by Gaia is expected to attain 2-3 fias for the stars with 
magnitude V < 10 mag and 10 fias for the stars of 
y = 15 mag. It is clear that not only the largest relativis- 
tic effects, but also many additional subtle effects should 
be taken into account to attain that accuracy. It is also 
quite clear that relativistic effects cannot be considered as 
small corrections to a Newtonian model as has been of- 
ten done earlier when the accuracy was not so high. The 
whole model should be formulated in a language compat- 
ible with general relativity. In such a relativistic frame- 
work many Newtonian concepts must be abandoned and 
the meaning of astrometric parameters such as position, 
parallax and proper motion of a star should be redefined. 



2. EXPERIMENTAL FOUNDATIONS OF GEN- 
ERAL RELATIVITY 

Einstein's general relativity is by no means the only pos- 
sible theory of gravity. However, it seems to be the sim- 
plest theory among the theories successfully passing all 
available observational tests. Let us briefly review the 
experimental foundations of general relativity. A detailed 
review of the modern experime ntal fo undations of gravi- 
tational physics can be found in Iwilil(l200lh . 



2.1. Einstein Equivalence Principle 

The basic principle of the theory is called Einstein Equiv- 
alence Principle. This principle consists of the following 
three parts: 



(1) The Weak Equivalence Principle stating that the 
masses on the both sides of the Newtonian gravitational 
law 

miner P = -G mgrav M r' / 

exactly coincide miner = mgrav for all bodies (actually, 
this is equivalent to the claim that G is a constant and 
its value is independent of the choice of the bodies with 
which we measure it). The Weak Equivalence Principle 
has been tested in many different experiments with a pre- 
cision of \dm\/m < 4 • lO^^'^. 

(2) Local Lorentz invariance stating that the outcome of 
any local non-gravitational experiment is independent of 
the velocity of the freely-falling test laboratory (refer- 
ence frame) where it is performed. This is equivalent to 
the principal postulate of special relativity theory which 
states that the light velocity in vacuum c is constant in 
any inertial reference system. This has been tested at a 
level of - 10^21 

(3) Local positional invariance which states that the out- 
come of any local non-gravitational experiment is inde- 
pendent of where and when in the universe it is per- 
formed. A part of local positional invariance can be tested 
by measuring of the gravitational red shift (e.g., of the 
clock frequency) 

Aiy/v=il + a) c-^ MJ, 

where a = in general relativity. A number of differ- 
ent experiments has proved that \a\ < 2 • 10^"*. Another 
part of local positional invariance (independence of "po- 
sition" in time) can be tested by looking for possible time- 
dependencies of fundamental (non-gravitational) con- 
stants. Different kinds of experimental data show tight 
constrains on possible time-dependence of the constants 
(e.g., the fine structure constant should be constant at 
a precision of — 10^^ over Hubble time of 13 billion 
years). 

2.2. Testing metric theories of gravity 

One can argue that if the Einstein Equivalence Principle is 
valid the gravity can be interpreted as an effect of curved 
spacetime. However, the Einstein Equivalence Principle 
does not necessarily imply general relativity. There exists 
a class of alternative theories of gravity compatible with 
that Principle. These theories are called metric theories 
of gravity. In order to test the principal observable ef- 
fects of metric theories of gravity a special scheme called 
Parametrized Post-Newtonian (PPN) formalism has been 
proposed (Will 1993). The scheme involves up to 10 nu- 
merical parameters which have different values in differ- 
ent theories and which can be fitted from observations. 
The most important PPN parameters are 7 and (3. Re- 
sults of data processing with a PPN reduction model in- 
colve a set of constraints on the PPN parameters. Alter- 
natively, the results can be interpreted as boundaries on 
possible deviations from general relativity. The most im- 
portant experimental results in this second interpretation 
are summarized in Table [2 

An additional test of general relativity is the search 
for possible time-dependence of Newtonian gravitational 
constant G. This can be done by looking for secular 
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Table 1. Most important post-Newtonian tests of general relativity 
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from general relativity 
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1 CI CI CllCC 


VLBI 


differential Shapiro delay 


±0.0003 


Eubanks et al. (1997) 


HIPPARCOS 


light deflection 


±0.003 


Froeschle et al. (1997) 


Viking radar ranging 


Shapiro delay 


±0.002 


^Reasenbera et al. (^222) 


Cassini radar ranging 


Shapiro delay 


±0.000023 


Bertotti et al. (2003) 


planetary observations 


perihelion advance 


±0.0002 


Pitieva (2001) 


Lunar laser ranging 


Nordtvedt effect 


±0.001 


Williams et al. (1996) 


Lunar laser ranging 


geodetic precession 


±0.001 


Williams et al. (1996) 



changes in semi-major axes of solar system planets (espe- 
cially. Mercury, Venus, Earth and Mars) as well as from 
pulsar timing of double pulsars. The most stringent es- 
timate here is \G/G\ < 10"^^ per year (jfitjeva 20Q^. 
General relativity predicts that G is time-independent. 
One more argument in favor of general relativity is the 
well-known indirect evidence for gravitational radiation 
in the double pulsar timing data. Gravitational radiation 
is strong-field regime phenomenon which is beyond the 
scope of the PPN formalism. 

All this shows that one can be "reasonably confident" 
about the correctness of general relativity, and that gen- 
eral relativity can be used as "standard" theory. Never- 
theless, it is still very important to test general relativity 
further. This will be discussed below in Section|6l 

3. RELAXmSTIC MODELING OF ASTRONOM- 
ICAL OBSERVATIONS 

Let us now outline general principles of relativistic mod- 
eling of astronomical observations. It is interesting that in 
spite of a deep conceptual difference between Newtonian 
physics and general relativity, the structure of the reduc- 
tion scheme changes, in principle, only in one point: light 
rays are no longer straight lines and should be carefully 
modeled. Figure |2l shows the four constituents of an as- 
tronomical observation in the relativistic framework. In 
curved spacetime there is no preferred coordinates where 
the laws of physics would have substantial simpler form 
than in other coordinates. Therefore, any reference sys- 
tem covering the spacetime region under study can be 
used. Instead of Newtonian inertial coordinates one has 
to choose some reference system in curved spacetime 
which is sketched symbolically on Figure |2] as a grid of 
curved coordinates. 

3.1. General scheme of relativistic modeling 

General scheme of relativistic modeling is presented on 
Figure |3l Starting from general theory of relativity, any 
other metric theory of gravity or the PPN formalism one 
should define at least one relativistic 4-dimensional refer- 
ence system covering the region of space-time where all 
the processes constituting particular kind of astronomi- 
cal observations are located. Each of four constituents 




Figure 2. Four parts of an astronomical event from the 
point of view of relativistic physics: 1) motion of the ob- 
served object; 2) motion of the observer; 3) trajectory 
of an electromagnetic signal from the observed object 
to the observer which represents a geodetic line (i.e. a 
"curved" line) in the chosen reference system; 4) the pro- 
cess of observation. The grid of curved coordinates in the 
background symbolizes the chosen relativistic reference 
system. 



of an astronomical observation should be modeled in the 
relativistic framework. The equations of motion of both 
the observed object and the observer relative to the cho- 
sen reference system should be derived and a method to 
solve these equations should be found. The equations 
of light propagation relative to the chosen reference sys- 
tem should be derived and a way to solve them should 
be found. The equations of motion of the object and 
the observer and the equations of light propagation en- 
able one to compute positions and velocities of the ob- 
ject, observer and the photon (light ray) with respect to 
that particular reference system at a given moment of the 
coordinate time, provided that the positions and veloci- 
ties at some initial epoch are known. However, the posi- 
tions and velocities calculated in this way obviously de- 
pend on the reference system, that is on the preferences 
of the person who writes down the equations. On the 
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Figure 3. General principles of relativistic modeling of 
astronomical observations (see text for further explana- 
tions). 



other hand, the results of observations cannot depend on 
the choice of the reference system. Therefore, it is clear 
that one more step of the modeling is needed: a rela- 
tivistic description of the process of observation. This 
part of the model allows one to compute a coordinate- 
independent theoretical prediction of observables start- 
ing from the coordinate-dependent quantities mentioned 
above. 

These four components can now be combined into rela- 
tivistic models of observables. The models give an ex- 
pression for relevant observables as a function of a set of 
parameters. These parameters can then be fitted to ob- 
servational data using some kind of parameter estimation 
scheme. The sets of certain estimated parameters appear- 
ing in the relativistic models of observables represent as- 
tronomical reference frames (see Section|5j- It is impor- 
tant to understand at this point that the relativistic mod- 
els contain some parameters which are defined only in 
the chosen reference system(s) and are thus coordinate- 
dependent. For example, position and velocity of the ob- 
served object are clearly coordinate-dependent. 



3.2. The Barycentric Celestial Reference System 

From the physical point of view any reference system 
covering the region of space-time under consideration 
can be used to describe physical phenomena within that 
region. In this sense we are free to choose the refer- 
ence system to be used to model the observations. How- 
ever, reference systems, in which mathematical descrip- 
tion of physical laws is in one sense or another simpler 
than in some other reference systems, are more conve- 
nient for practical calculations. Therefore, one can use 
the freedom to choose the reference system to make the 
parametrization as convenient and reasonable as possible. 

Two Working Groups on relativity in astrometry, celestial 
mechanics and metrology established 1997 by the Inter- 



national Astronomical Union (lAU) and Bureau Interna- 
tional des Poids and Mesure (BIPM) have come to the 
conclusion that the most convenient relativistic reference 
system for the applications in astrometry, solar system 
dynamics, and time keeping an d dissemination is defined 
by the following metric tensor (ISoffel et al .120031) : 



ffoo 
9tj 



2w 

72" 



2w^ 

^4 



+ Oic-% 



rW 



0(c- 



l + ^w]+ 0{c 



(1) 



with the post-Newtonian potentials w and defined by 

w{t,x)^G [ d'x'^^^ 



2c2 dt^ 



/.3.V(,.')I-.'I,(2) 
(3) 



a and cr* being related to the components of the energy- 
momentum tensor T"'^: 



1 



r 



lOO 



1 



T 



(4) 



The origin of spatial coordinates of this reference sys- 
tem is chosen to coincide with the barycenter of the so- 
lar system. The reference system defined in this way is 
called Barycentric Celestial Reference System (BCRS). 
The BCRS has been explicitly recommended by the lAU 
for the modeling of high a ccuracy astr onomical obs erva- 
tions (lAU 2001; Rickman 2001;'Soffel ef flZ."2003'). For 
the moment the BCRS is a post-Newtonian reference sys- 
tem with higher order terms (post-post-Newtonian terms, 
etc.) neglected in the metric tensor Q. The reason for 
that is that the post-Newtonian approximation is suffi- 
cient to model any observations in forseeable future (in- 
cluding microarcsecond astrometry as long as the obser- 
vations are made further than about one degree from the 
Sun). Post-post-Newtonian terms can be added to the 
metric tensor as soon as they are necessary for some ap- 
plications. The word "celestial" in the name of BCRS is 
used to underline that the BCRS do not rotate with the 
Earth and that remote sources do not move relative to the 
BCRS in some averaged sense. The second reference sys- 
tem deined by the same lAU resolutions (Rickman 2001) 
is the Geocentric Celestial Reference System (GCRS). 
This reference system is only marginally important for 
Gaia (mostly for modeling of orbit trac king data and re - 
lating the Gaia onboard clock to TAI (Klioner 2003a")) 
and will not be discussed here. The PPN version of the 
BCRS valid for certain class of metric theories of gravity 
can be found in lKlioner & Soffeil(l200(l andlWiill Il993l) . 

The BCRS will be also used for the modeling of Gaia 
observations. This is a reference system underlying the 
resulting Gaia catalogue (see Section|5lbelow). The co- 
ordinate time of the BCRS is called Barycentric Coordi- 
nate Time (TCB). The TCB will be used to parametrize 
the Gaia catalogue. 
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3.3. Motion of the objects and the observer 



can be written as follows 



Typically, for objects situated in the solar system (aster- 
oids, planets, space vehicles) the equations of motion are 
ordinary differential equations of second order and nu- 
merical integration with suitable initial or boundary con- 
ditions can be used to solve them. For objects outside 
of the solar system one use often simple models like uni- 
form and rectilinear motion in space or more complicated 
ones, e.g., for binary stars. In any case one should under- 
stand that in the relativistic framework all these ad hoc 
models give positions and velocities of observed objects 
in the chosen relativistic reference system. 

The principal relativistic effects in the translational mo- 
tion of bodies in the solar system (including Gaia satel- 
lite, asteroids, etc.) are contained in the so-called 
Einstein-Infeld-Hoffmann (EIH) equations of motion of 
N gravitating bodies, whose gravitational fields can be 
described by their masses Ma only: 



GMb 



XA - Xb 



\xa - Xb\ 



\Fpn{Mb,xb,xb) + 0(0" 



(5) 



The Newtonian part of these equations (shown explicitly 
above) follows from the term of order in goo- The rel- 
ativistic terms require all other terms in the BCRS met- 
ric tensor specified above. Various parts of these equa- 
tions represent: (1) relativistic perihelion advance (^43 
"/cty for Mercury, ~10 "/cty for Icarus, etc.); (2) geode- 
tic precession (~2 "/cty for Lunar orbit); (3) various pe- 
riodic relativistic effects (important mostly for LLR and 
binary pulsar timing observations). Further effects not 
contained in the EIH equations are the effects due to ro- 
tation of the bodies (Lense-Thirring or gravitomagnetic 
effects) and those due to non-sphericity of the gravitating 
bodies. These additional effects are marginal for the cur- 
rent accuracy of LLR and SLR, but negligible for Gaia. 
In case of Gaia satellite one should use a slightly simpli- 
fied version of the EIH equations since the influence of 
the mass of the satellite on the motion of other gravitat- 
ing bodies can be neglected. 

The BCRS metric tensor allows one also to derive the 
equations of rotational motion of an extended body. 
These equations will not be discussed here, since they 
are not important for Gaia. 



3.4. Light propagation 

In any metric theory of gravity the equations of light 
propagation coincide with the equations of geodetic lines 
in the chosen reference system. The latter are ordinary 
differential equations of second order These equations 
could also be solved by numerical integrations, but nor- 
mally one prefer to use some approximate analytical so- 
lutions. Only in some special (normally, highly symmet- 
rical) cases like Schwarzschild metric exact analytical so- 
lutions are known. Anyway, an appropriate way to solve 
the equations of light propagation should be found. 

The structure of the BCRS equations of light propagation 



x{t) = Xq + ccr (t - to) + c 



(6) 



where Xq and cr are the parameters of Newtonian straight 
line, SpAT are the post-Newtonian terms, and Si^^p^ are 
the additional effects induced by the motion of gravitat- 
ing matter (i.e., by translational and rotational motion of 
gravitating bodies). The terms of order of c^^ in both 
goo and gij are required to derive SpAr(i), and the terms 
in goi are needed for Si^^p^. The next order effects, 
the so-called post-post-Newtonian effects, would require 
terms of order of c^'' in both goo and gij (the c^'* terms 
in gij are not in the current definition of the BCRS metric 
tensor). The principle observable effects in the light prop- 
agation are (1) the gravitational light deflection (amount- 
ing to 1.75" for a light ray grazing the Sun) and (2) the 
gravitational signal retardation (the Shapiro effect; this 
effect amounts to ^ 240 fis for the radar ranging of Venus 
in upper conjunction). 



3.5. Conversion to observables: proper direction 

As mentioned above the conversion the coordinate- 
dependent quantities into coordinate-independent ob- 
servables is an important part of relativistic modeling. 
From the mathematically point of view the coordinate- 
independent quantities are scalars. Special mathematical 
techniques are known to perform the suitable conversion 
in each particular case. One of the most important ap- 
plication of this conversion procedure is a conversion of 
the coordinate direction n into the source into the cor- 
responding observable direction s. The observable di- 
rection is often called "proper direction" in gravitational 
physics. Proper direction is a direction relative to the 
proper reference frame of the observer (see Section [S] 
about the difference of the concept of "reference frame" 
in astronomy and gravitational physics). A proper ref- 
erence frame is a mathematical model of an ideal clock 
and three orthogonal rigid rods which the observer uses 
to measure time intervals, distances and directions in his 
vicinity. In special theory of relativity the proper refer- 
ence frame of an observer is related to some inertial ref- 
erence system by a Lorentz transformation. It is there- 
fore, sufficient to use Lorentz transformations to convert 
n into s. The parameter of the Lorentz transformation in 
this case coincides with the velocity of the observer rela- 
tive to the chosen reference system. In general relativity 
it is also sufficient to use Lorentz transformations, but the 
parameter 1/ of the transformations should be related to 
the BCRS velocity of the observer as 



l + —w{t,Xo) ] +0(c-'*), 



(7) 



where Xo and Xq are the BCRS positions and velocity of 
the observer, respectively. A detailed discussion of this 
conversion and comparison of different approaches can 
be found in Klioner (2003b). The relativistic terms in 
(Q are derived from the c^^ terms in goo and g^ of the 
BCRS metric tensor The difference between n and s 
can be called relativistic aberration. The difference be- 
tween the Newtonian aberration and the relativistic one 
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Figure 4. Five principal vectors used in the model {see 
text for explanations). 

may amount to several milliarcsecond for Gaia observa- 
tions. 

3.6. Conversion to observables: proper time 

Another important case is the conversion of intervals of 
the coordinate time t into the corresponding intervals of 
the proper time r of the observer The general form of 
this conversion reads 

dT /dt=l + ApN + c-^ AppN + 0(c-5), (8) 

where ApM and App^ are the post-Newtonian and post- 
post-Newtonian terms, respectively. Explicit form of 
these two functions depends on the metric tensor: in or- 
der to compute for Apjv the c^^ terms in goo are needed, 
while the c^"* terms in goo^ the c^^ terms in goi. and the 
c^^ ones in are required to compute ^ppjv- Typically 
in the Solar system and in particular for Gaia onboard 
clocks Ic-Mpjvl - 10-^ and \c-'^ AppM \ ^ lO^i^. 

4. RELATIVITY FOR GAIA 

Now, having all these theoretical tools one can formulate 
the relativistic model for Gaia. The relativistic model for 
Gaia is well documented (Klioner 2003a), so that we just 
outline the overall structure of the model here. 

4.1. Structure of tlie standard relativistic model 

The model consists essentially in subsequent transforma- 
tions between 5 following vectors (Figure|3l: 

a) s is the unit observed direction (the word "unit" means 
here and below that the formally Euclidean scalar product 
s ■ s = is equal to unity). 



b) n is the unit vector tangential to the light ray at the 
moment of observation, 

c) <T is the unit vector tangential to the light ray at t = 

— oo, 

d) k is the unit coordinate vector from the source to the 
observer, 

e) I is the unit vector from the barycenter of the Solar 
system to the source. 

Note that the last four vectors should be interpreted as 
sets of three numbers characterizing the position of the 
source with respect to the BCRS. Vector s represents 
components of the observed direction relative to the local 
proper reference system of the satellite. All these vec- 
tors would change their numerical values if some other 
relativistic reference system is used instead of the BCRS. 
The model consists then in a sequence of transformations 
between these vectors as shown on Figure [S] The phys- 
ical meaning of each transformation can be summarized 
as follows (the numbering here coincides with the num- 
bering on Figure|5|i: 

(1) aberration (effects vanishing together with the 
barycentric velocity of the observer): this step converts 
the observed direction to the source s into the unit BCRS 
coordinate velocity of the light ray n at the point of ob- 
servation; 

(2) gravitational light deflection for the source at infinity: 
this step converts n into the unit direction of propagation 
(T of the light ray infinitely far from the solar system at 

t —>■ —(x; 

(3) coupUng of finite distance to the source and the grav- 
itational light deflection in the gravitational field of the 
solar system: this step converts cr into a unit BCRS coor- 
dinate direction fe going from the source to the observer; 

(4) parallax: this step converts k into a unit BCRS direc- 
tion I going from the barycenter of the Solar system to 
the source; 

(5) proper motion, etc: this step provides a reasonable 
parametrization of the time dependence of I (and, possi- 
bly, of the parallax tt) caused by the motion of the source 
relative to the barycenter of the Solar system; 

(6) orbit determination process. 

These transformations ha ve alr eady been discussed in 
full detail (ISioner 2003at iKlion er & Peit^ 2003; Klione] 
I2003tf). Let us only mention the following. The most 
complicated part of the model is the light deflection 
model where the effects of ( 1 ) monopole fields of all ma- 
jor solar system bodies, (2) quadrupole fields of the gi- 
ant planets, and (3) gravitomagnetic fields due to trans- 
lational motion of all major bodies should be taken into 
account in order to attain the accuracy of 1 /las. More- 
over, each body with a mean density p and radius R > 
(p/1 g/cm"^)"^/^ X 650 km produces a light deflection 
of at least 1 //as. Therefore, a few tens of minor bod- 
ies (mainly, satellites of the giant planets) should also be 
taken into account in certain rare cases (K lioner 2003a) . 
The parametrization of time dependence of I in the rela- 
tivistic framework looks exactly the same as in the New- 
tonian case. The only difference is that all vectors and 
parameters here (parallax, proper motion, etc.) are coor- 
dinate quantities defined in the BCRS. 
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remote sources: 



(1) (2) (3) , (4) , (5) s . s , . 

s < > n < > <T < > k < > I, TV < > lyLQ), 7r(toJ) ^^['to), 



1 1 • (1) (2.3) , (6) , 

solar system objects: s < — > n ^—>- k < — > orbit 

Figure 5. Transformation sequences (see text for explanations). 



4.2. Implementation of the model 

An ANSI C code has been written to imple- 
ment the relativistic model in its full co mplexity 
(iKhoner & Bla nkenburg 2003; Klioner 2003b). The 
model has been implemented in two modes: predictor 
mode and corrector mode. Predictor mode implements 
the standard way of astrometric reductions when the ob- 
served direction to the source is predicted starting from 
some a priori catalogue parameters (coordinates, proper 
motion, parallax, etc.) of that source. The catalogue is 
supposed to be improved later by fitting the parameters 
to the whole set of data. Corrector mode implements the 
reductions in the opposite direction, that is, the momen- 
tary barycentric direction to the source is restored from 
the observed direction as good as possible. The code does 
not contain any attempt to restore parallaxes and proper 
motions, or orbits of the sources: the model can only be 
applied separately for each individual observation. 

Both modes are implemented both for solar system ob- 
jects and for remote sources situated outside of the solar 
system. The principal difference between remote sources 
and solar system sources lies in the treatment of the 
light propagation: an initial value solution of the corre- 
sponding differential equations is used in the former case, 
while in the latter case two point boundary value problem 
should be solved. Although analytical approximations 
are used in both cases, a rather time-consuming numer- 
ical inversion process is used for solar system sources in 
the predictor mode in order to attain the goal accuracy 
of 1 /ias. The corresponding refinement of the model 
aime d at direct analytical solution for this case is under- 
way llKlioner & Blankenburgl2004 . 

For remote sources the corrector and predictor modes 
being implemented independently of each other 
must give exactly the opposite transformations. 
This was used to massively test the implementation 
(Klioner & Blankenburs 2003). The situation is different 
for solar system objects. Here for the corrector mode, it 
is statistically better to calculate the gravitational light 
deflection as if the body were a remote source, even if it 
is known a priori that the source is a solar system body 
(but it is not known how far the body is, otherwise at 
least a preliminary orbit is known and one should better 
use the predictor mode). 

The implementation was conceived to be as flexible as 
possible. Internal parameters allow one to select the type 
of arithmetic to be used (to test possible numerical in- 
stabilities), to change easily any of the physical, math- 
ematical and astronomical constants used in the model 
(including switching between several available planetary 
ephemerides), to switch on and off each individual effect. 
Both predictor and corrector mode routines have a goal 
accuracy parameter, which is used together with some a 
priori criteria to decide which effects should be computed 
in each particular case. The latter feature allows one to 



speed up the calculations substantially if a lower accu- 
racy is sufficient (e.g. the source brightness information 
can be used to meet the Gaia observational accuracy for 
fainter objects). 

The implementation has been used in massive numeri- 
cal tests aimed at identifying possible inconsistencies or 
numerical instabilities as well as points of critical numer- 
ical performance. The implementation was used also to 
test anot her simpli fied model implementation used in the 
GDAAS JAn glada -Escude et al. 2003|. The implementa- 
tion of the full model will be further supported and opti- 
mized. 



4.3. Beyond the standard relativistic model 

The model described above is constructed under assump- 
tion that the solar system is isolated. This means that any 
influence of gravitational fields generated outside of the 
solar system are ignored in the model. For the majority 
of the sources the external field can indeed be fully ne- 
glected, but there are a number of cases when the exter- 
nal gravitational fields produce observable effects. Sev- 
eral authors ha ve discus s ed thes e additional effects in de- 
tail (s ee, e.g., iKlioneJ ^2003a^ and Kopeikin & Gwini] 
(12000) '). Let us briefly list here the main effects of this 
kind: 

(1) Gravitational light deflection caused by the masses 
situated outside of the solar system: (a) weak microlens- 
ing on the stars of the Galaxy (Belokurov & Evans 2002), 
(b) lensing on gravitational waves (both primordial ones 
and those from compact sources), (c) lensing of the com- 
panions of edge-on binary systems. 

(2) Cosmological effects. 

(3) More complicated models for the motions of observed 
objects in the BCRS are necessary for the case of binary 
stars, etc. 

Note that all these effects can be easily taken into ac- 
count by a simple additive extension of the standard 
model since at the required accuracy the external grav- 
itational fields can be linearly superimposed on the solar 
system gravitational field. The only exception could be 
the effects of cosmological backgroun d, but a prelimi- 
nary study by iKhoner & Soffell (l2004 >) shows that even 
here the coupling of the local solar system fields and the 
external ones can be neglected. 



5. GAIA REFERENCE FRAME 

It is important to remember that all astrometric parame- 
ters of sources obtained from Gaia observations will be 
defined in the BCRS coordinates: positions, proper mo- 
tions, parallaxes, radial velocities, orbits of minor plan- 
ets, binaries, etc. All these parameters will represent the 
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Gaia reference frame, which is a materialization of the 
BCRS. The Gaia reference frame is, so to say, a model of 
the universe in the BCRS. Thus, the goal of astrometry in 
the relativistic framework is not to find "the" barycentric 
inertial reference system, which is unique in Newtonian 
formulation, but to find a materiahzation of some chosen 
relativistic reference system. 

Let us note here that the meaning of words "reference 
system" and "reference frame" in relativistic astronomy 
is different from the meaning normally used in gravita- 
tional physics. Reference system is a purely mathemat- 
ical construction (a chart) giving "names" to space-time 
events. A reference frame is, in contrast, some material- 
ization (realization) of a reference system. In astronomy 
the materialization is normally given in a form of a cata- 
logue (or ephemeris) containing positions of some celes- 
tial objects relative to the selected reference system. Any 
astronomical reference frame (a catalogue, an ephemeris, 
etc) is defined only through the reference system(s) used 
to construct physical models of observations. 



6. GAIA FOR RELATIVITY 

Using general relativity for the standard reduction model 
does not mean that Gaia data should not be used to test 
general relativity itself. On the contrary, testing relativ- 
ity is one of the exciting goals of the mission. Gaia will 
certainly deliver an estimate of the PPN parameter 7, ap- 
pearing mainly in the magnitude of the light deflection 
effects, with an unprecedented accuracy of '--^ 5 x 10^^. 
However, it is by no means the only way Gaia will im- 
prove our knowledge of gravitational physics. Gravita- 
tional light deflection could be tested in a much more 
profound way with the Gaia data Gaia. It will be cer- 
tainly possible to look for terms with totally different de- 
pendence on the angular distance to the deflecting body. 
In this sense one should be able to get first experimen- 
tal estimates of higher-order effects. Also the PPN pa- 
rameter (3, appearing in the equations of motion of so- 
lar system bo dies, will be determined w ith an accuracy 
of - 10-4 (Hestrofer & Berthier"2004) which IS com- 
parable with the current accuracy from the planetary fits 
(Pitieva 2001). A simultaneous fit with the planetary 
data could further improve the accuracy. Special data 
processing of observations close to Jupiter and Saturn 
(which should be dropped from the global absolute so- 
lution since the positions of those planets are not known 
with an accuracy necessary to predict the light deflec- 
tion at the level of 1 /ias) will allow to test subtle rel- 
ativistic effects caused by translational motion of the 
planets and by their quadrupole gravitational fields (see 
ICrosta & Mignard (2004) for a preliminary study of the 
second of these possibilities). A number of cosmolog- 
ical tests (upper estimates of parallaxes of quasars, ap- 
parent proper motions of quasars, possible traces of low- 
frequency gravitational waves in those proper motions, 
direct measurement of the acceleration of the solar sys- 
tem relative to quasars, etc.) will also have a big impact 
on our knowledge. A lot of work is still necessary to un- 
derstand how to use the full potential of the huge amount 
of observation data Gaia will deliver to us in the most 
efficient and useful way. 
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